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We have computed the quark anti-quark expectation value in finite volume at two loop in chiral perturbation 
theory and compare it with a formula obtained in analogy to the Liischer formula for pion mass in finite volume. 
We observe that due to the small finite size correction at two loop it is not possible to obtain conclusions on 
the accuracy of the extended Liischer formula. 



1. Introduction 

Quantum chromo dynamics (QCD) is a strongly 
coupled theory at low energy. Therefore it is difficult 
to study the physics below the QCD scale of about 
1 GeV. Lattice QCD on the other hand is an attempt 
to tackle this problem by taking the Lagrangian and 
calculating (parts of) the generating functional nu- 
merically. Although these computations with fairly 
small quark masses are now possible, at the same time 
the finite volume effects arc also becoming important. 
In this work we employ chiral perturbation theory 
(ChPT) in order to study these effects analytically. 
Chiral perturbation theory, in its modern form, was 
first proposed by Weinberg Q and extended in 0, . 
ChPT is an effective field theory to describe the strong 
interactions at low energy. To do the computation in 
finite volume as suggested in Q one imposes a periodic 
boundary condition on fields which result in the mo- 
mentum quantization and consequently modification 
on the quantum corrections. Since the applicability of 
ChPT is limited to small momenta, in finite volume 
this leads to 



F n L > 1 . 



(1) 



which holds for L larger than about 2 fm. F„ is the 
pion decay constant and L is the linear size of the 
boundary. In addition we stay in a limit where the 
Compton wavelength of the pion is smaller than the 
lattice size which corresponds to 



mlFlV » 1 



(2) 



This is the so-called p-regime. This condition ensures 
that perturbative treatment will work. Expectation 
value of the quark condensate per se will be achieved 
by varying the generating functional with respect to 
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the scalar external field as follows 
1 S 

< qT a q > = _(_i_)Z e// [« o ]| ia=0 . (3) 

Where T a stands for the Gell-Mann matrices. We 
then carry out the (qq) at two loop in finite volume 
within the framework of chiral perturbation theory. 
This quantity has already been evaluated at one loop 
order in Another way for the computation of the 
finite volume effect is the Liischer approach. In this 
method one obtains the leading finite size effects for 
the pion mass to all order in perturbation theory from 
the scattering amplitude [EQ- It has later been ex- 
tended to evaluate the volume dependence of the pion 
decay constant 8]. As the second goal we shall then 
obtain a formula for the the quark condensate volume 
dependence by following the same line of reasoning 
as in the Liischer work. This work was published in 
Ref. @. 



2. A Liischer formula for the vacuum 
condensate 

We were inspired by the work of Liischer, where 
he showed how the finite size effect on the masses of 
spin-less particles are related to the forward elastic 
scattering amplitude in infinite volume. This occurs 
when particles are enclosed in a lattice box and as 
considered in this letter with temporal direction of 
the space-time extended to infinity. To see explicitly 
where these effects stem from, one looks at the 
two-point correlation function with fields subjected 
to a periodic boundary condition. The correlator 
takes on the form 



(4) 
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(a) (b) 

Figure 1: Two potentially relevant graphs to the (qq). 
® indicates the insertion of (qq) and the solid line is the 
meson propagator 

where the inverse propagator reads 

G L (f,m%)=p 2 -m a -^ L { P 2 ) (5) 

and Yjl(jj 2 ) denotes the sum of the one-particle irre- 
ducible diagrams to be evaluated in momentum space 
with discretized momenta and p = 2 jp- , with ft a three 
dimensional vector of integers. In practice the desired 
quantity to compute is the change in £(p 2 ) due to 
the finite volume correction and the difficult task in 
Liischer's work was to show that this consists of all 
self-energy diagrams with just one propagator in each 
diagram allowed to wrap around the whole position 
space and therefore experiencing the boundary. Using 
the Poisson summation formula the modified propaga- 
tor can be obtained in terms of the propagator in the 
infinite volume accompanied by an exponential factor 
as follows 

G L (q) =J2 e ~ i9iniLiG °° 

n 

where in three dimensions the length squared for a 
given ft reads k — ft 2 — n\ + n| + n 2 ,- For ft = 
it gives back the propagator in infinite volume. The 
crucial point comes about in Euclidean space where 
<7o — * ~ «<?o and, as shown below, in our result for (qq) 
the exponential factor falls off rapidly for large ft 2 . 
Liischer in his analysis kept only ft 2 = 1 but summa- 
tion over all values shall be carried out numerically in 
our formula as suggested in |10| . 

Out of three self-energy graphs shown in Ref. @, 
only two are potentially relevant to the (qq) as de- 
picted in Fig. ^ The diagram (a) has however no 
contribution when we apply ChPT to the case of par- 
ity even operators. Extending our task to the evalua- 
tion of a general operator in the the diagram (b) we 



arrive at the following relation 
<0)y-(0)oo= / ShJ2 e ~ lLq - nG ^\°\^- ( ? ) 

J [ ' n=iO 

Integration over momentum can be split into paral- 
lel and transverse components with respect to the ft, 
followed by the continuation qn — ► qu — is. For s — > oo 
the integration over this path is negligible. Unlike the 
mass case it turns out that the computation of the di- 
agram (b) for the (qq) is simpler mainly because the 
matrix element (</>\0\(fi) representing the amputated 
vertex function has no external momenta in it and in 
fact it should be evaluated at zero momentum trans- 
fer. After isolating the poles and picking out the right 
one we, complete the integration over which finally 
leads to the following formula 



(O v )-(0 00 ) = -(cf>\0\ ( f>) 




As the last step we identify the remaining integral as 
the modified Bessel function K\ 

(O) v -{0) 00 = -{cj>\0\<l>) E 16^2 /77M ^ l(C(fc)) 

(9) 

where Q(k) = y/(k)mL and x(k) is the multiplicity,!. e. 
x(k) is the number of times the value k = ft 2 appears 
in the sum overt ft. Deriving this formula also has the 
advantage that one can estimate the size of the sigma 
term by evaluating the volume dependence of the the 
quark anti-quark expectation value. 

3. The vacuum expectation value at two 
loop 

The calculation of the quark condensate at two loop 
was done in 01 . We go back to these computation 
while taking into account the finite size effects. As 
mentioned before, this effect comes about as a modi- 
fication in the propagator. In Fig. [21 we show all the 
Feynman diagrams that contribute up to p e . In infi- 
nite volume the loop diagrams we deal with, contain 
the functions 

B <"' 2) = (4/ <, V^W <10) 

These integrals are solved in the dimensional regu- 
larization scheme with e = (4 — d)/2 . 
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Figure 2: Feynman diagrams up to p 6 for (qq). 
vertices o, and gi stand for the insertion of (qq) at p 2 , 
p 2 and p 6 respectively. • is the p 2 and x the p 4 vertices. 
Solid lines indicates the meson propagator 
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A(m 2 ) 



B(m 2 ) = 



16tT 2 



7 B + log(±n) + 1 



+A{m 2 ) + eA t {m 2 ) + o{( 2 ) 
- 7e + logiiTT) + 1 



(11) 



16tt 2 



+B(m 2 ) + eB e (m 2 ) + o(e 2 ). (12) 



At finite volume integration over momenta should 
be replaced by a summation over discretized mo- 
menta. It is also important to note that the pole 
structure does not change its form, since the singu- 
larity arises when momenta in the loop become very 
large and therefore the pole remains unchanged de- 
spite the replacement of integration with summation. 
The final result for the integrals at finite volume be- 



A(m 2 



m 2 ( m 2 
16^ l0g ( -3 



1 



16tH 



B(m 2 



1 



16tt 2 
1 

16tt 2 



fe=i, 



log 



c(ky 



/J,* J 167T 2 

E x(k)2K (C(k)) (14) 



k— l,oo 



with Qik) — VkmL. It is evident from equations 
(O and lfH|) that the simple relation A(m 2 ) = 
m 2 (l/167r 2 + B(m 2 )) is no longer valid in finite vol- 
ume computations. After putting all the diagrams 
together and applying the dimensional regularization 
scheme with renormalization parameter fi we ensure 
that the divergent parts cancel and for finite part we 
use the expressions in (|13|l and (|14fl instead of those 
in t| 1 1 1> and i|12|l . The parts containing A e and B e 
cancel in the final result for the vacuum expectation 



Figure 3: The ratio R u is shown for three different pion 
masses 91 . 



value. We denote the lowest order masses for pions, 
kaons and etas as Xw> Xk and Xn respectively. They 
are given in terms of the strange quark mass, m s and 
the average of up quark and down quark masses, m 
by 



2B rh , xk = B (m s + rh) 
X-q = 2B {2m s + m)/3 . 
We also define the following quantities 



(15) 



- -B F 2 



1 



-B Q F 2 1 



(«) 

V 4 

F 

v\ 
~Fh 



F 2 



F 2 



The analytical expressions for 
evaluated in ChPT and can be found in 



,(«) 



(16) 



Vq are 



. Numerical 

results are shown in Fig.[3]and Fig.^J where the ratio 



R q is defined as 



Ra 



(qq)v - (qq)c 



(17) 



The C\, the low ener 
from diagrams at p 6 



y coupling constants (LECs) 
cancel out in the numerator 
and we have ignored them in the denominator but the 
L\ will be considered as inputs and we use the values 
from fit 10 of 12]. As expected the volume size effects 
are pronounced for small L and small pion mass. But 
the full NNLO contribution to the finite volume effect 
is much smaller than the total NNLO numerical cor- 
rection. In turn this total NNLO correction is smaller 
than the one presented in [llj . This is due to the fact 
that in Ref. [ll| the physical value of pion decay con- 
stant Fj[ and also physical masses for mesons are put 
in the formula. 
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Figure 4: The ration R s is shown for three different pion 
masses increasing from top to bottom [?|. 



4. Conclusions 

We compared the quark condensate in finite volume 
obtained in two different ways. The method referred 
to as Lrischer's approach yields the leading term of 
this quantity to all orders in perturbation theory while 
only one propagator experiences the finite volume. 
One the other hand, we try to do this computation 
up to p e rather directly by letting the modification 
of all propagators be present at the same time. The 
comparison of the direct calculation at one-loop with 
the extended formula is trivial but at two-loop level, 
however, this is not the case. This allows in principle 
a check of the accuracy of the extended Liischer ap- 
proach. Unfortunately, because of the very small size 
of the finite volume effect at p 6 , no conclusion can be 
drawn concerning the accuracy level of the extended 
Liischer formula from our work. 
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